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Abstract. We show that the function P{E) derived from the density of states of a 
constant heat capacity reservoir coupled to some system of interest is not identical to 
the physically measurable (transitive) temperature. There are, however, connections 
between the two quantities as well as with the Tsallis parameter q. We exemplify these 
connections using the one-dimensional Ising model in the "dynamical ensemble" . 
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There is an ever growing list of systems that are shown — or beheved — to be inadequately 
described by standard Boltzmann-Gibbs statistical mechanics. All of these systems 
in some way involve one or more of the following common features: (i) long-range 
interactions, (ii) long-time memory, (iii) ( mult i-) fractal space-time. In order to tackle 
these systems mathematically, Tsallis [IJ introduced a one-parameter generalization 
of Shannon's information entropy and, along the lines of Jaynes j2], established the 
framework of generalized statistical mechanics (GSM) jSlllj- 

GSM has been successfully applied to a wealth of problems — for a regularly updated list 
see jS]. There is, however, to this day some confusion about the appropriate definition 
of the thermodynamic temperature. Furthermore, the physical interpretation of the 
nonextensivity parameter q is — in most cases — still an open question. 
This paper is organized as follows: In section 1 we shall give a short introduction to 
Tsallis entropy, generalized mean values and canonical probability distributions. In 
section 2 we shall introduce the "physical temperature" and the possible transformation 
to an extensive entropy. We shall show that this transformation is of little consequence 
for practical calculations. In section 3 we summarize how the "Tsallis factor" can be 
obtained from a microcanonical basis in the case of a "constant heat capacity" reservoir, 
and section 4 is devoted to an application of Tsallis statistics to the nearest-neighbour 
Ising model in contact with an ideal gas. Our conclusions will be given in section 5. 
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1. Tsallis statistics 

The generalized entropy introduced by Tsallis (and previously in different contexts by 
Havrda, Charvat and Daroczy) reads 

where i = 1, . . . , f2 labels the possible microstates of the system under consideration, {pi\ 
are the microstate probabilities, g is a real parameter and is a positive real constant 
that approaches Boltzmann's constant k-Q as Sq approaches the Shannon entropy in the 
limit g ^ 1: 

Sq'^ Si = S = -k^i^Pilnpi. (2) 

i=l 

In the following fcx will be set to unity. 

The parameter q determines the degree of nonextensivity (nonadditivity) of Sq in the 
following sense: If two statistically independent systems A and B with probabilities 
|pf I and |pf |, respectively, are combined to form the system AVJ B, the composite 
entropy Sq^^ satisfies the relation 

Sf^ = + 5f + (1 - g) S^S^ . (3) 

Thus, the entropy is superextensive (superadditive) for g < 1, extensive (additive) for 
g = 1 and subextensive (subadditive) for g > 1. 

Sq has been shown to be positive for all g, to exhibit definite curvature for g < 
and g > (it is constant for g = 0), and to attain its maximum value in the case of 
equiprobability. 

The connection to thermodynamics has been established by generalizing the canonical 
ensemble to the case g 7^ 1, i.e. maximizing the entropy under the constraints 
of (i) normalization of the probability distribution and (ii) knowledge of the energy 
expectation value. 

The normalized g-expectation value jlj of some observable O is defined as 

{{Oh = %fV ' 

where {Oj} are the microstate values of O. 
Maximization of the functional 



leads to the generalized canonical probability distribution 



1 



3TMP 



l-q 



(z = l,...,n), (6) 
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where TMP stands for Tsallis-Mendes-Plastino to distinguish the parameter (3 from 
other choices (introduced below), Zq is the generahzed partition function, 



i=l 



TMP 



2^j=i Vj 



1-9 



(7) 



and [. . .]+ symbohzes the Tsalhs cut-off condition, i.e. [. . .]+ = [. . .]©(. . .), where Q{x) 
is Heaviside's unit step function. 

With this formahsm, the Legendre transform structure of thermodjTiamics is preserved 
IHj, in particular 

In order to solve the manifestly self-referential equations (jHl) for the probabilities pi, 
Tsallis et al. introduced the so-called transformation |^: Using an auxiliary 

parameter (3', equation (jHl) can be rewritten in the form 

p.^ (9) 



e;l, ii - (1 - <i)P'E,]x 



where 



Martinez et al. J\ showed that with equations ^ for the microstate probabilities a 
maximum of Sq cannot be guaranteed since the concomitant Hessian is not diagonal. 
To overcome this problem, they suggest a reformulation of the internal energy constraint, 
thus maximizing the functional 



TMP 



F' [fe}] = Sq-a(j:p.-l]- {E, 

\i=l ) i=l 



m,) ■ (11) 



Here OLM stands for "Optimal Lagrange Multiplier". This procedure leads to the 
entropy maximizing probabilities 

1 



Vi 



^OLM 



(l-q)P^^^(E,-{{E))q) , (12) 



which can also be written in the form Q with 

^OLM 

^' = l + (l-g)/?°^^((i?)), ■ ^^^^ 

In this case the Legendre structure is altered, and equation (jH)) reads 

Various other choices of the internal energy constraint have been made, among those the 
reformulation of the variational problem in terms of the so-called "escort probabilities" 
Pi = Pi/ J2j jHIHl- Most interesting, however, is a work by Suyari P, who showed that. 
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if one wants to interpret Sq as the average of some nonextensive information measure 
Iq and postulates that Iq and 5*^ satisfy the same nonextensivity relation, one not only 
has to use normalized g-expectation values, but also the "normalized" g-entropy 

^norm ^'1=1 fi fl5) 

The same "normalization" of the Tsallis entropy has previously been suggested by 
Rajagopal and Abe '^^ grounds of form invariance considerations of the generalized 
KuUback-Leibler relative entropy. 

As will be shown later in this paper, the actual choice of the form of the energy constraint 
is insubstantial as far as one (i) uses normalized g-expectation values and (ii) is only 
interested in the behaviour of the system as a function of the "physical temperature" 
to be introduced in the following section. 



2. Physical temperature 

Consider an isolated system G composed of two subsystems A and B in "weak" thermal 
contact. "Weak" contact means that the interaction energy between the two systems 
shall be negligible, thus 

Eg = Ea + Eb = const. (16) 

If A and B do not have the same temperature from the start, after some (possibly very 
long) time, thermal equilibrium will be reached. This equilibrium state is characterized 
by a maximum of the composite entropy under the constraint of constant total energy. 
Since the entropy Sq^^ satisfies equation Q, one gets [TT] 



+ 



d{{EAh 



B 



which, because of 6{{EA))q = —S{{EB))q, leads to the equilibrium parameter ("physical 
temperature" ) 



— = Z^^phys 



1 ( \ _ 1 / 



1 + (1 - q)S^ \d{{EA))qJ 1 + (1 - q)S^ \d{{EB))qJ 

This result is sometimes referred to as the "generalized zeroth law of thermodynamics" . 
For the two cases from the previous section this means that 







TMP 



/3phys = -^r-^ = . (19) 



The first identity in (jl9p has also been obtained by Kalyana Rama |12 on grounds of 
the equivalence between canonical and microcanonical ensemble. 

The necessity of the transformation ()18|1 has also been shown by Vives et al. ^3] from 
quite a different point of view: Using equation Q and scaling of the entropy, they 
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established a generalization of the Gibbs-Duhem relation, showing that the Lagrange 
parameters (in the TMP formalism) are not intensive quantities. By introducing the 
above transformation of the Lagrange parameters together with a transformation of the 
entropy, 



the standard Gibbs-Duhem relation can be reobtained. 

The same transformations have been deduced by Johal ^3] from the existence of thermal 
equilibrium of two systems in weak contact and the existence of a Legendre transform 
structure. The new entropy Sq is an extensive quantity, i.e. S^^^ = + S^, and by 
inserting equation (|T)) one can see that Sq is the Renyi entropy: 



Sq = Sn^ , ^=^/^ . (21) 



1-q 

As has been shown by Lenzi et al. extremization of Sn with the constraints 

of normalization and normalized g-expectation value of the energy leads to the OLM 
distribution (fT^ with 

dS 

P = /^Phys = . (22) 

The complete Legendre transform structure is preserved jTHl UH] • 

Now, it is interesting to note that for actual calculations the probabilities (jH)) or 
(|12p must be determined either iteratively or by means of the (3-f3' transformation. 
The latter procedure entails calculation of the probabilities in terms of the auxiliary 
parameter /?', computation of the relevant averages and subsequent determination of the 
associated temperature jHITTj. If one is interested in the behaviour of the system under 
consideration with the "physical temperature", one gets for all formulations (TMP, 
OLM, normalized entropy, Renyi entropy) the same relation: 

^■^'^^ ^ 1 - (1 - q)P'{{E))q ■ ^^^^ 
However, if one wishes to address the subject of thermodynamic stability, some subtleties 
have to be considered. The generalized free energy that is minimized by the OLM 
distribution ()12|) is associated by Legendre transformation not with the Tsallis entropy, 
but with the Renyi entropy [Tn|. It reads 

i^^""" = {{Eh - l^Sn = {{E))q - In ^ . (24) 

Pphys Pphys 



Defining 



lnZOL^=lnZOL^-/5ph,s((i?))., (25) 



thermodynamic relations recover their familiar forms ^3] 



Pphys 
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^-^ phys f^-^ phys 

Since thermodynamic stability is commonly associated with positivity of the specific 
heat, one gathers from equation ()28j) that F^^^ is required to be a concave function 
of the physical temperature. As will be illustrated in section 4, due to the /3phys~ 
/?' transformation, this criterion is not always fulfilled, and some kind of Maxwell 
construction must be applied. Also, for reasons of concavity of the Renyi entropy S^, 
the parameter q is restricted to the interval < g < 1. 

3. Microcanonical derivation of the Tsallis factor 

As in the previous section we consider an isolated system G composed of two weakly 
interacting systems S ("System") and B ("Bath"). Once again, the energies satisfy the 
relation 

Eg = Es + Eb = const. (29) 

The principle of equal a-priori probabilities demands that, in a situation where the total 
energy is Eq and the system energy is Es (the bath energy is accordingly Eg — Es) , all 
compatible states are equally probable, thus 

f-'*^"'" ikiiEiP, • 

where Qx{E) is the number of states (the structure function) of system X with energy 
E, and VLg{Eg) is given by 

nG{EG) = / ns{Es)nB{EG - Es) dEs . (31) 

Jo 

From the definition 

(3{E) = ^ In^BiE) (32) 
the Boltzmann factor 

Ps{Es) oc exp{-pEs) (33) 
can be obtained, provided one of the following conditions holds: 

• The temperature of the bath is exactly constant ( "constant T" derivation) or 

• The energy of the bath is large compared to the system energy: Es <^ Eb ("small 
Es" derivation) 
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The latter condition, however, does not unambiguously lead to Boltzmann-Gibbs 
statistics JH]- Indeed, every distribution of the form 



Ps{Es) cx exp2(-/5E5) (34) 

can be obtained, where expQ(x) = (1 + QxY^^ is a Q-deformed exponential that 
converges to the ordinary exponential in the limit Q 0. Of this, the "Tsallis factor" 
is but a special case (for Q = 1 — q). 

By an extension of the "small i?^" derivation to the case where Es is not necessarily 
small (i.e. there is no truncation of some series expansion involved), Almeida fOl EH! 
has derived the above distribution (jH^ under the assumption that the heat capacity of 
the reservoir is exactly constant: 

AEb 1 d / 1 \ ^ 

^' = ^=Q= ^ JeS [Wb) ) = ^ = ^^'^ 

This includes the true heat reservoir limit, where the temperature is exactly constant, 
irrespective of the amount of heat gained or lost by the bath. In this case the heat 
capacity of the bath is infinite, thus Q = 0, which is the Boltzmann case. Note that 
with (|35|) the parameter Q and therefore the Tsallis parameter q is linked to a physical 
property of the heat bath and thus gets a physical interpretation. So q is restricted to 
g < 1, if one does not allow for negative heat capacities. 

Starting from equation together with (jS21), Tatsuaki [THj established the following 
differential equation for the structure function Qb{Eb) of the bath: 

dlnQniEB) = . ^ n° ^Eb , (36) 

where f3o = (3{0) is an integration constant. Integrating this equation, one finds the 
solution 

nBiEB) = ^BiO) expQiPoEB) . (37) 



From this and equation ()32j) one gets for the P function of the constant heat capacity 
environment 

Using Eb = Eg — Es, one obtains 

^BiEc - Es) = ^BiO) expQ [(3oEg] exp^ [-P{Eg)Es] , (39) 

and, inserting this into equation (jHUI) and using (jHT|) . one finally arrives at 

p.iEs) = ^siEs)e.p^[-PiEG)Es] 

J,^^nsiEs)expQ[-/3iEG)Es]dEs' ^ ' 



with 
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Comparison with equation Q shows that the two probabihty distributions are the same 
with 

/3(i?G) = /3', (42) 
and from equation one gets that the (inverse) physical temperature is given by 

where = {{Eq - Es))q = Eq - {{Es))g has been used. 

For the sake of completeness, let us mention that already in 1994 Plastino and Plastino 
PT] showed that for a system in contact with a finite heat bath with structure function 
n{E) oc E°' the resulting equilibrium distribution is of Tsallis form with q = {a — l)/a. 

4. The Ising model in the "dynamical ensemble" 

As an application, let us now consider a one-dimensional Ising model (S) of spins with 
periodic boundary conditions coupled to an ideal gas {B) with M degrees of freedom 

m- 

The Ising model is given by the Hamiltonian 

N 

H = J2il- a,a,+,) , (44) 
1=1 

where (Tj is a spin variable with values ±1, corresponding to the spin pointing "up" or 
"down" . The total energy of the composite system is given by 

Eg = Es + Eb = Ne + Mk = const. , (45) 

where e = Es/N is the energy per spin of the Ising system, and k = Eb/M is the kinetic 

energy per degree of freedom of the gas. 

With the known structure function of an ideal gas. 



M~2 



nBiEB) = CEs' , (46) 

C being some constant and 

Eb = Eg — Es , (47) 

we get for the (unnormalized) probability of finding the composite system G in a state 
where S has the energy Es: 

M-2 

Ps{Es) cxns{Es){EG-Es)— . (48) 
Since the ideal gas is a constant heat capacity environment, 

we expect the equilibrium distribution of the Ising system to be of Tsallis form with 

2 , , 
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From equation (j32j) follows 



ms) ^M^^Ui-l). (51) 



2Eb 2k \ MJ 

One sees that in the limit of an infinite heat bath M — > cxd the statement of the 
equipartition theorem is recovered. In that limit, (3{Eb) is the (inverse) physical 
temperature of the bath. 

To be able to calculate the energy of the Ising system as a function of the physical 
temperature by means of Monte Carlo simulations, it is necessary to determine the 
normalized g-expectation value of the system energy {{Es))q at given total energy Eq- 
This is accomplished using the generalized Metropolis sampling algorithm by Andricioaei 
and Straub [23 123- Weighting the microstate energies E^^ with p| (-E's'*) , the transition 
probabilites from state i to j, related by single spin flips, become 



P{i j) = min 



M-4- 
iJ) \ — 



(52) 



where equations (pSj) and (jSUj) have been used. From equations and we obtain 
the common physical temperature of the bath and the Ising system: 

^Phy. = (Eg - {{Es)),) . (53) 

Figure 1 shows the energy per spin of an = 128 Ising system as a function of the 
physical temperature for different sizes of the gas reservoir. The data points were 
obtained by means of Monte Carlo simulations, while the lines represent exact results, 
calculated from equations ^ and using the known structure function of the Ising 
model and the appropriate q values given by (jSHI)- The curve for g = 1 is also given for 
comparison. 

Besides the fact that the simulation and the exact results are in good agreement, one 
finds that the energy is a multi-valued function of the physical temperature. Thus, we 
take a look at the generalized free energy defined by ()24|) and discover a behaviour that 
is similar to what was described by Lima et al. in for the TMP case. Figure 2 shows 
the free energy (a) and the internal energy (b) for the case q = ji {M = 24). The solid 
lines represent _pOLM {{e))g calculated by means of the Pptys^P' transformation, while 
the broken curves were calculated from the iteratively determined probabilities 
The /?phys^/3' transformation produces a closed loop in the free energy curve, whereas 
the iterative method reveals metastable branches. The behaviour depends on whether 
one starts from low or high temperatures. By always choosing the branch with the low- 
est free energy the metastable states are eliminated, leading to a discontinuous internal 
energy as indicated in figure 2 (b). 
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Figure 1. Specific energy of an A'' = 128 Ising system for bath sizes M = 18 and 
M = 34 as a function of the physical temperature. Symbols represent simulation 
results, lines are exact results obtained from the g-canonical distribution for the 
appropriate q values. 

5. Conclusions 

We have shown that the function P{Eb), defined in analogy to the Boltzmann-Gibbs 
{q — 1) microcanonical ensemble as the logarithmic derivative of the structure function 
of a constant heat capacity reservoir B coupled to some system S, does not describe 
the common physical temperature of the two systems. Instead, I3{Eb) at Eb = Eg-, 
where Eq is the constant total energy of the system+bath, corresponds to the auxiliary 
parameter (3' introduced in Tsallis canonical ensemble theory. From this the physical 
temperature can be obtained by a /3phys^/3' transformation (which is the same for all 
formulations of the generalized canonical ensemble) , provided one knows the normalized 
^'-expectation value of the bath energy Eb at constant Eq- 

Using the one-dimensional Ising model in contact with a finite ideal gas reservoir as 
an example, we have shown that this transformation can be performed adopting the 
generalized Metropolis Monte Carlo sampling algorithm of Andricioaei and Straub. 
The curves of the free and internal energy thus obtained, show the same "reentrant 
behaviour" already known from the TMP canonical ensemble. Choosing always the 
states with the lowest free energy, single-valued curves can be obtained. The curves 
of the internal energy become discontinuous, while the free energies have discontinuous 
first derivatives. 
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Figure 2. Free energy (a) and internal energy (b) of an = 128 Ising system for 
q = Yj. The solid lines are obtained via the Pphys~P' transformation, the broken 
curves are calculated iteratively. Latter curves have been slightly displaced for better 
visualization. 



In summary we can say that a physical system that is in weak contact with a constant 
heat capacity reservoir is described by a Tsalhs canonical distribution with an index 
q that is determined by the nature of the bath (including the case of an infinite heat 
capacity, in which the Boltzmann-Gibbs distribution is recovered). The factor f3' in the 
g-exponential is only the inverse temperature of the system if g = 1. 



Acknowledgments 

The author wishes to thank Prof. Dr. G. Miinster for the support and fruitful discussions 
during the preparation of this paper. 



[1] Tsalhs C 1988 J. Stat. Phys. 52 479 

[2] Jaynes E T 1957 Phys. Rev. 106 620 

[3] Curado E M F and Tsalhs C 1991 J. Phys. A 24 L69 

[4] Tsalhs C, Mendes R S and Plastino A R 1998 Physica A 261 534 

[5] http:/ /tsallis. cat. cbpf.br/biblio. html 

[6] Yamano T 2000 Eur. Phys. J. B 18 103 

[7] Martinez S, Nicolas F, Pennini F and Plastino A 2000 Physica A 286 489 



Physical temperature and the meaning of the q parameter in Tsallis statistics 



Abe S 2000 Phys. Lett. A 275 250 
Suyari H 2002 Phys. Rev. E 65 066118 
Rajagopal A K and Abe S 1999 Phys. Rev. Lett. 83 1711 

Abe S, Martinez S, Pennini F and Plastino A 2000 Preprm^ |cond-mat/0011012 
Kalyana Rama S 2000 Phys. Lett. A 276 103 
Vives E and Planes A 2002 Phys. Rev. Lett. 88 020601 
Johal R S 2002 Preprint cond-niat/0207268 

Lenzi E K, Mendes R S and da Silva L R 2000 Physica A 280 337 
Abe S 2001 Physica A 300 417 

Lima A R and Penna T J P 1999 Phys. Lett. A 256 221 



Wada Tatsuaki 2002 Preprint cond-mat/0208205 



Almeida M P 2001 Physica A 300 424 

Almeida M P, Potiguar F Q and Costa U M S 2002 Preprint [cond-mat /0206243 
Plastino A R and Plastino A 1994 Phys. Lett. A 193 140 
Gerling R W and Hiiller A 1993 Z. Phys. B 90 207 
Andricioaei I and Straub J E 1996 Phys. Rev. E 53 R3055 
Andricioaei I and Straub J E 1997 Physica A 247 553 



